Nonlocal phenomenology for anisotropic MHD turbulence 
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A non-local cascade model for anisotropic MHD turbulence in the presence of a guiding magnetic 
field is proposed. The model takes into account that (a) energy cascades in an anisotropic manner 
and as a result a different estimate for the cascade rate in the direction parallel and perpendicular to 
the guiding field is made, (b) the interactions that result in the cascade are between different scales. 
Eddies with wave numbers fcii and k± interact with eddies with wave numbers q\\ , q± such that a 
resonance condition between the wave numbers q\\,q± and kn,k± holds. As a consequence energy 
from the eddy with wave numbers ku and k± cascades due to interactions with eddies located in the 

resonant manifold whose wavenumbers are determined by: q« ~ e 1 / 3 ^ 3 / B, q± = k± and energy 
will cascade along the lines k\\ ~ C + k 2 / 3 e 1 ^ 3 /Bo- For a uniform energy injection rate in the parallel 
direction the resulting energy spectrum is E(k\\ ,k±) — e 2 ^ 3 k^ 1 k ± 5 ^ 3 . For a general forcing however 
the model suggests a non-universal behavior. The connections with previous models, numerical 
simulations and weak turbulence theory are discussed. 

PACS numbers: 



I. INTRODUCTION 

Magnetic fields are met very often in astrophysics; in- 
terstellar medium, accretion discs, the interior of stars 
and planets. In most of these cases the magnetic fields are 
strong enough to play a dynamical role in the evolution of 
the involved astronomical objects [2(j. Furthermore, the 
fluid and magnetic Reynolds numbers are large enough 
so that a large number of scales are exited and coupled 
together, making it very difficult to calculate the evolu- 
tion of these systems even with the power of present day 
computers. As a result a turbulence theory that models 
the behavior of the small unresolved scales is in need. 
The simplest set of equations that describes the evolu- 
tion of the flow and the magnetic field when the two are 
coupled together are the magneto-hydro-dynamic (MHD) 
equations that in the Elsasser formulation are written as 



d t z ± = ±B ■ Wz ± - z T ■ Vz* - VP + z/V 



2 z± 



(1) 



where z ± = u ± b with u the velocity, b the magnetic 
field, v the molecular viscosity assumed here equal to the 
magnetic diffusivity r\ — v, B is a uniform magnetic field, 
and incompressibility V • z = has been assumed. 

For zero viscosity the above equations conserve three 
quadratic invariants the magnetic Helicity (that we are 
not going to be concerned with in the present work) and 
the two energies E ± — J(z ± ) 2 dx 3 . The question then 
arises in the limit of infinite Reynolds number is there a 
physical process under which the two energies cascade to 
sufficiently small scales so that they can be dissipated? 

For hydrodynamic turbulence a description of such a 
process exists and was given by Kolmogorov 13] (K41). 
In his phenomenological description the energy zf at a 
scale / interacts with similar size eddies and cascades in 
a timescale l/zi. As a result in a statistically steady state 
the energy cascades in a scale independent way at a rate 



nomenological description of the energy cascade in hydro- 
dynamic turbulence there have been attempts to derive 
similar results for MHD flows. However, non-trivial diffi- 
culties arise when a mean magnetic field is present. First 
the MHD equations are no longer isotropic resulting in 
an anisotropic energy flux and energy spectrum. Sim- 
ple dimensional arguments can not be used to estimate 
the degree of anisotropy that is a dimensionless quantity. 
Second the MHD equations are no longer scale invari- 
ant, as a result simple power law behavior of the energy 
spectrum is expected only in the small or large B limit 
that scale similarity is recovered. Finally, it is not clear 
that interactions of similar size eddies (local interactions) 
dominate the cascade. Different size eddies could play an 
important in cascading the energy. 

The first model for MHD turbulence was proposed by 
Iroshnikov pj] and by Kraichnan [14] (IK). The IK- 
model assumes isotropy and that the time scale of the 
interactions of two wave packets of size I is given by the 
Alfven-time scale t a ~ l/B. The energy cascade due to 
a single collision is given by Az 2 ~ (z 3 /1)t a . The num- 
ber of random collisions that would be required then to 
cascade the energy is going to be N ~ (z 2 /Az 2 ) 2 . As 
a result the energy will cascade in a rate e ~ E/Nt a ~ 
z 4 /(Bl) and therefore z t ~ (eB) 1 ^! 1 ^. The resulting 1- 
D energy spectrum is then given by E(k) ~ (eB) 1 / 2 ^ 3 / 2 . 
The assumption of isotropy however has been criticized 
in the literature and anisotropic models have been pro- 
posed for the energy spectrum. Goldreich & Sridhar 
[lfj ] (GS) proposed that in strong turbulence the cas- 
cade happens for eddies such that the Alfven time scale 
t, ~ Bkn is of the same order with the non-linear time 
scale t nl ~ zk± (so called critical balance relation), 
where fcii and k± are the parallel and perpendicular to the 
mean magnetic field wavenumbers respectively. Repeat- 
ing the Kolmogorov arguments then one ends up with 



// that leads to the prediction z t ~ l 1/3 or in terms the energy spectrum E(k\\ , k±) — fc ± 5/3 with the parallel 
of the 1-D energy spectrum E(k) ~ fc -5 / 3 . Since the phe- and perpendicular wave numbers following the relation 
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fcii ~ . A generalization of this result was proposed 
by Galtier et al. Q where the ratio of the two time scales 
t a /t nl was kept fixed but not necessarily of order one, in 
an attempt to model MHD turbulence both in the weak 
and the strong limit. Bhattacharjee & Ng Q (BN) re- 
peated the IK-modcl arguments replacing the nonlinear 
time scale by t nl ~ l±/zi and the Alfven time scale by 
t a ~ l\\/B. Further assuming that the cascade is hap- 
pening only in the k± direction obtained the energy spec- 
trum E(k) ~ (eBY^kJ 2 ^ 1 . Finally Zhou, Matthaeus 
& Dmitruk [2l| (ZMD) suggested using as time scale the 
one given by the inverse average of the Alfven and non- 
liner time scale r _1 = (r. "p 1 + r _1 to obtain a smooth 
transition from the K41 to the IK and the anisotropic BN 
result depending on the amplitude of the guiding field. 

Although this large variety of models exists the agree- 
ment with observations II ll and with the results of nu- 
merical simulations 0, i,Hl, M, G3, 01 is only partially 
satisfactory and seems to be case dependent. Further- 
more, all these models assume locality of interactions 
(i.e. only similar size eddies interact and only one length 
scale is needed in each direction in the phenomenological 
description). Locality of interactions however has been 
shown to be in question by both theoretical arguments 
and analysis of data in numerical simulations even in the 
isotropic case [l], H, Iflj . Furthermore the rigorous result 
for weak turbulence @, [j| has shown that only modes 
in the resonant manifold = are responsible for the 
energy cascade. It seams reasonable therefore that non- 
locality is an essential ingredient of MHD turbulence that 
needs to be taken in to account in a model. In addition 
we expect that the energy will not cascade isotropically 
so not only the amplitude of the energy cascade rate is of 
importance but also the direction. With these two points 
in mind (anisotropy and non-locality) we try to construct 
a non-local model, for the energy cascade. 

To begin with the new model let us consider a MHD 
flow in a statistically steady turbulent state forced at 
large scales, in the presence of mean magnetic field 
B. We will denote the two 2-D energy spectra as 
E + (k±, fcii) and E~(k±,k\t) where the total energy is 
given by E^ = J E^dkxdkn. For simplicity assume 
E~ (fc_L, fen) ~ E + (k±, k\\) (i.e. negligible cross helicity) 
and drop the ± indexes leaving the case E~ <C E + to be 
investigated in the future. To shorten the notation we 
will write Ek = E(k±, fen). The index k denotes that E k 
depends on the wavenumbers fc|| and k±. 

Let us now consider two eddies of different scales z£ 
and z~ interacting. Let us assume that the z^ eddy 
has wave numbers ~ k±,ku and the z~ eddy has wave 
numbers ~ Qx,q\\< Here we will focus on the cascade of 
the energy of the z£ eddy, the cascade of the z~ can be 
obtained by changing the indexes k, q and ±. From the 
form of the non-linear term we expect by dimensional 
analysis that the rate of energy cascade of the z^ eddy 



will be 

£(k) ~{zt) 2 {z-)q 

~[k ± k\\E k ][q m E q ]V 2 q (2) 

and similar for the energy cascade rate of the z~ eddy. 
Note that in such an interaction the energy will not cas- 
cade isotropically but it will depend on the value of q. In 
a interaction of the two eddies the energy of the z£ eddy 
will move from the wavenumber k to the wave number 
k + q. If <7|| -C q± then most of the cascade will be in 
the q± direction. As a result we need to define separately 
the rate energy cascades to larger k±: £±(k) and the rate 
energy cascades to larger h\ : £\\ (k) as: 

£ ± ^[k ± k ll E k ][q ± q ll E q } 1 / 2 q ± (3) 

£\\ ~ [k x k\\E k ][q ± q ] \E q ] 1 %. (4) 

Note that in writing the equations above we have not 
taken in to account possible scale dependent correlations 
between the two fields that could reduce the energy cas- 
cade. Such an effect has been taken into account by 
Boldyrev based on the GS model and could be in- 
corporated in the present model. However we will not 
make such an attempt here since we want to present the 
model in its simplest form. Equations [3] and [4] express 
the rate energy cascades in the absence of a mean mag- 
netic field and are valid only when \q\ < \k\ because small 
eddies although they have a stronger shear rate z q q decor- 
rclatc making them less effective in cascading the energy. 
However in the presence of guiding field not all wave 
numbers q are as effective in cascading the energy E k . 
Because the two eddies z k ,z q travel in opposite direc- 
tions the time they will interact will be the Alfven time 
t a ~ [(7||£?o] _1 and the time needed to cascade the energy 
will be t nl ~ (z^l)- 1 ~ Hqly/qxq^Eq}- 1 . Therefore, 
from all the available wavenumbers only the wave num- 
bers with t a > t nl will be effective in cascading the 
energy. This restriction leads to: 



1\\B < q±yq± q\\E q (5) 

where we used the approximation \q\ ~ q± as a first order 
approximation of \q\ for large B. This relation looks very 
similar to the critical balance relation of the GS model. 
However in this case the relation [5] gives the wave num- 
bers that the eddy z£ will interact with and does not 
restrict the location of z^ in spectral space. We are go- 
ing to refer to the set of wavenumbers that satisfy the 
relation above as the resonant manifold and we are going 
to use this relation as an equality because from the al- 
lowed wavenumbers qu the ones closer to ku will be more 
effective in cascading the energy. Here as in the equations 
[3]and|4l we assume that q\\ < fcy. Finally since the mean 
magnetic field does not directly effect the _L direction we 
assume that similar size k± and q± are the most effective 
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FIG. 1: An illustration of the non local model. A eddy at 
(ku,k±) interacts (non-locally) with an eddy at the resonant 
manifold (gn (k±), k±) cascading the energy in the direction of 
the arrow. 



to cascade the energy q± ~ kx (i.e. locality in the k± 
direction). So equation [5] is written as 



k\E q 

R2 
D 



(6) 



We are now ready to impose the constant energy flux 
condition that would lead to a stationery spectrum. Be- 
cause the cascade is anisotropic, constant energy flux now 
reads: 



d kl £± + 0fc..£|| 







(7) 



The equations I3I4I6|7I form the basic equations of our 
model. It is worth noting, that in this model the cascade 
of energy decreases by the introduction of the guiding 
field not because the individual interactions weakened 
but because the number of modes that are able to cascade 
the energy decreases due to the resonance condition [51 A 
sketch of the mechanisms involved in the model is shown 
in figure [1] 

First let us concider the weak turbulence limit that 
is obtained in the limit B — > oo. For large B based 
on equation [5] the resonant manifold becomes very thin 
q\\/k± <C 1. Furthermore because £\\/£x ~ q\\/kx <C 1 
we can neglect the cascade in the parallel direction. If 
also E k is non-singular at fcii = we have that E q ~ 
E(kx,0). Substituting qn from [5] in [3] and imposing the 
constant flux condition [7] we obtain: 



£± ~kik ll E k E{k ± ,0)/B = e(k ll ). 



(8) 



Note that the spectrum E(k±,k\\) depends on the en- 
ergy of the resonant manifold E(k±,0) just like the weak 
turbulence result and unlike what the BN local theory 
for weak turbulence predicts. If the energy spectrum for 
E(kx,kt\) and for the resonant manifold E(k±,0) scale 
like k\ and k™ respectively then we end up with the 
weak turbulence prediction 



y_tl 



n = -4 



(9) 



Assuming that the two spectra are smooth around fc|| = 
(just like the weak turbulence theory needs to assume) 
we obtain: 



E" 



(10) 



As we decrease the value of B we need to take in to 
account that the energy cascade in the || direction is non 
zero. In this case energy does not cascade in the _L di- 
rection but cascades along the lines that are tangent to 
the direction off and satisfy dk\\/dkx = £^/£\\ = q\\/k± 
or 



k± Q\\(k' x : 



dk' ± + C. 



(11) 



Let A be the length along such a curve; then we can move 
to a new coordinate system given by (A, C). In this new 
coordinate system the constant flux relation to first order 
in qn reads 



or 



^ = A.{[k ± k {l E k }[q xqi] E q n\ 



[k 1 _k\\E k \[k^E q ] 1 ' 2 k x = e{C) 







(12) 



(13) 



where only terms up to order q» are kept. Letting — » qu 
we obtain the equation for the resonant manifold. 

[k x q\\Eg][kxq\\Eg]V 2 k x = e(0) = e . (14) 

Substituting qn from [6] and solving for E q we obtain 

E q = E(q l{ (k x ), k X ) = el /3 Bkl 7/3 . (15) 
The equations [5] and [TT] then gives us 



„2/3 1/3, 



| - „_l *4 /a /B and fc|, = ^k 2 / 3 el /3 /B + C. (16) 

Returning to the equation for the energy energy [13] we 
get 



E k 



<C)e- 1/3 

kTh 



(17) 



where C is given by equation 1111 and the predicted spec- 
tra 1 171181 are valid in the range q\\ < k\\ < oo. For smaller 
values of k\\ the condition /c|| < q» that we initially as- 
sumed is not satisfied. The energy of the modes inside 
the resonant manifold is given by [15] and no singularity 
at fey = exists. 

In the special case that e(C) is a constant e(C) = eo 
that corresponds to a uniform injection rate per unit of 
wavenumber (fen) at the large scales k± — > the spectrum 
reduces to 



2/3 



E k 



1 - 5 /3, 
k ± fc„ 



(18) 
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but in general the spectrum will depend on the way en- 
ergy is injected in the system. The non-universality that 
the model suggests is due to the fact that we assumed 
that the energy cascades in a deterministic way only 
along the lines in the (k\\,k±) plane given by [IT] In re- 
ality energy will not cascade strictly along the lines [TTJ 
but there is going to be some exchange of energy be- 
tween lines that could bring the energy spectrum in the 
form [TBI However, if and how fast a universal spectrum 
can be obtained in MHD is not an easy question to an- 
swer. This question is related to the return to isotropy 
of an anisotropicaly forced flow in hydrodynamic turbu- 
lence, that is still an open question. If indeed in MHD 
turbulence in the presence of a guiding field there is a 
universal spectrum this is expected to happen at smaller 
scales than in hydrodynamic turbulence because nonlin- 
ear interactions are weaker. 

So far we concerned ourselves with only large values 
of B. In principle we could extend our results to smaller 
values of B with out making some of the approximations 
used to arrive at the results 1 171 181 However, such pro- 
cedure leads to more complex equations that prevent us 
from deriving the energy spectrum in a compact form 
and we do not make such an attempt at present. 

It is worth emphasizing the similarities the current 
model has with GS model. Both models emphasize the 
role of the manifold fen ~ k\ Cq 3 /B, obtained by the 
resonance condition [5] in out model or the critical bal- 
ance condition in the GS model. However in this model 
the cascade is not restricted in this manifold, but instead 
all modes in the (feii,fej_) plane cascade due to non-local 
interactions with the modes in this manifold. 

An other point we need to emphasize is that taking 
B — ► oo does not reduce the predicted spectra in [T7] to 
the weak turbulence limit [TUJ This is due to the two 



different limiting procedures followed. In the first case 
feq fTT)]) first the limit B — > oo was taken and then the 1-D 
flux was determined, while in the second case [T7] we first 
obtained the 2-D energy flux and then the limit B — > oo 
was taken. Note however that the condition [9] is satisfied 
in both cases and the spectrum is also smooth at ku = qii 

— 7/3 

because although the resonant manifold scales like k ± 
while the rest of the spectrum scales like k^^k^ 1 the 
resonant manifold widens as k± increases. It is possible 
as we discuss bellow that in different (numerical) setups 
either of the two limiting procedures can be valid and 
different spectra could be obtained. 

In numerical simulations a finite discrete number of 
modes is kept. Based on this model the cascade rate 
is reduced in the presence of a mean magnetic field not 
because the individual interactions themselves arc weak- 
ened but because the number of modes that that interact 
effectively is reduced. If the modes in a numerical simu- 
lation with the smallest non-zero wavenumber k\ = 2n/L 
(where L is the box height) is larger than the resonant 
manifold ( k\ > q\\ ) then if B is further increased the 
scaling with B of the energy dissipation rate will be lost 
and the spectrum exponents could change, since the num- 
ber of modes in the resonant manifold already have taken 
their minimum value (i.e. the number of modes that have 
feii = 0). Furthermore a difference in the energy spectrum 
exponents can be expected in numerical simulations if the 
modes inside the resonant manifold are not forced. The 
sensitivity of the model to the way the system is forced 
could in part explain the disagreement in the measured 
spectrum exponents. 
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